Unconventional Quantum Hall Effect and Tunable Spin Hall Effect in M0S2 Trilayers 
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We analyze the Landau level (LL) structure and spin Hall effect in a M0S2 trilayer. Due to 
orbital asymmetry, the low-energy Dirac fermions become heavily massive and the LL energies grow 
linearly with rather than with V^. Spin-orbital couplings break spin and valley degenerate LL's 
into two time reversal invariant groups, with LL crossing effects present in the valence bands. We 
find a field-dependent unconventional Hall plateau sequence v — ■ ■ ■ —2M — 6, — 2M — 4, — 2M — 2, 
— 2M — 1, • • • , —5, —3, —1, 0, 2, 4 • • • . In a p-n junction, spin-resolved fractionally quantized 
conductance appears in two-terminal measurements with a controllable spin-polarized current that 
can be probed at the interface. We also show the tunability of zero-field spin Hall conductivity. 



PACS numbers: 73.43.-f, 71.70.Di, 72.25.Mk, 75.76. +j 

Introduction. — Successful isolation of a single molyb- 
denum disulfide (M0S2) trilayer presents a new platform 
to explore interesting two dimensional (2D) electronic 
physics [IHl]- The honeycomb lattice structure of a M0S2 
trilayer (when viewed from the top) and its low-energy 
Dirac physics are reminiscent of graphene [5]. Indeed, the 
M0S2 trilayer exhibits advantages over graphene in sev- 
eral areas of intense recent interest, specifically, whether 
this 2D material has an energy gap [6]-[8] and whether it 
has substantial spin-orbit couplings (SOC) [9UTT]. Unfor- 
tunately, the answers for graphene to date are still not 
satisfactory, even though tremendous efforts have been 
made to improve the possibilities [121 US] • 

When the layered compound M0S2 is thinned down to 
a single trilayer, it departs from an indirect gap material 
to a direct gap material [2l O [T4HT7] . Additionally, some 
density functional theory (DFT) calculations have shown 
that there exists large SOC in M0S2 pH2Q]. The exis- 
tence of a large energy gap and strong SOC may place 
this newly discovered 2D platform ahead of graphene in 
the race for the next generation of semiconductors. More 
recently a M0S2 transistor with room-temperature mo- 
bility about 200cm^/(V-s) has already appeared [4l. 

In this Letter, for the first time, we analyze how SOC 
infiuences the Landau level (LL) spectrum of massive 
Dirac fermions and how to increase the spin Hall ef- 
fect in a M0S2 trilayer. We find that the Hall plateau is 
field-dependent and follows an unconventional sequence 
that has an even-odd-even transition. The LL energies 
grow linearly with 5, and the spin and valley degenerate 
LL's are broken into two time reversal invariant groups. 
The broken symmetry in the valence band also gives 
rise to LL crossing effects, leading to pronounced peaks 
in the measurement of longitudinal magnetoresistance. 
We further investigate the case of a p-n junction, where 
spin-resolved fractionally quantized conductance appears 
in two- terminal measurements, with a controllable spin- 
polarized current that can be probed at the interface by 



scanning tunneling microscope (STM). We also explic- 
itly show how to tune and increase the zero-field spin 
Hall conductivity by reducing the inversion asymmetry. 
None of these exotic features is able to be observed in 
graphene, even with the help of electron-electron inter- 
actions. 

Continuum theory. — We start from a description of 
the low-energy model of an isolated M0S2 trilayer, which 
applies generally to other group- VI dichalcogenides with 
the same crystal structure. The top and bottom S layers 
and the middle Mo layer are parallel triangular lattices. 
Because of their ABA relative stacking order, the top 
view of this trilayer forms a honeycomb lattice with S 
and Mo atoms at A and B sites, respectively. Near the 
Brillouin zone inequivalent corners K and K', the con- 
duction and valence band states are approximately from 
|</>c) = \dz^) and = {\d^2_y2) ^ iTz\d^y))/V^ or- 

bitals, respectively. This effective two-band model has 
been suggested by DFT calculations [20^ and supported 
by optical experiments [2TH24] . To linear order in p, the 
effective k-p Hamiltonian in the above basis reads 

n = v{pxTz(Jx + Py(^y) + ^(7^ - Xt^SzCTz + At^S;^ , (1) 

where the Pauli matrices cr operate on the space of the 
and d^id orbitals, = ±1 labels the K and K' val- 
leys, and Sz = ±1 denotes the electron spin t and |. The 
Fermi velocity v is given by V = at/h 0.53 x 10^ m/s, 
where t is the effective hopping between the two d- 
orbitals of Mo mediated by the p-orbitals of S. As 
anticipated, the inversion asymmetry [25 between d^^ 
and d ± id orbitals gives rise to the Aa^ mass term 
which pins the ground state to a quantum valley Hall 
(QVH) insulator [26l|27]. Mo atoms provide strong in- 
trinsic SOC ~ TzSzCFz O [27] that adjusts the energy 
gaps to 2(A — A) for TzSz = 1 bands and to 2(A + A) 
for TzSz = —1 bands. Note that this SOC perturba- 
tion preserves inversion (V = TxCTx) [25] and time re- 
versal (T = iTxSyK) symmetries. As a combined ef- 
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feet of breaking inversion symmetry and strong SOC, the 
T invariant term Ar^s^ breaks the particle- hole symme- 
try by oppositely shifting the r^s^ = ±1 bands. Using 
A = 830 meV and A = 37.5 meV extracted from DFT cal- 
culations [20 , Fig. [l] plots the band structure of a M0S2 
trilayer which exhibits two features that substantially dif- 
fer from graphene, i.e., (i) a large QVH band gap and (ii) 
a lifted degeneracy between TzSz = ±1 bands. We ad- 
dress that while the conduction band bottoms are lined 
up for all flavors, the valence band tops have a signifi- 
cant shift in energy between TzSz = ±1 bands. These 
symmetry breaking features in the spin- valley space can 
be further verified by the flavor-dependent energy disper- 
sions 
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E± XtzSz ± 



(A - XtzSz 



(2) 



where ± stands for the conduction and valence bands. 

Broken symmetry LL's and LL crossing effects. — In 
the presence of a uniform magnetic field perpendicular to 
a M0S2 trilayer, the 2D kinetic momentum p in Eq. ([T]) 
is replaced by tt = p + eA/c. We now derive the LL 
spectrum using the Landau gauge A = (O^Bx). The 
operators tt = tt^ -\- iTVy coincide with the lowering oper- 
ators, satisfying Trc/)^ = —i{h/iB)V^<l^n-i and nt^o = 0. 
Here £3 = V^/{eB) = 25.6/^/B\f] nm is the magnetic 
length and is the nth LL eigenstate of an ordinary 
2DEG. This model is approximately valid when Hv/Ib 
is smaller than the band width ~ 500 meV. To focus on 
the influences from SOC and inversion asymmetry on the 
LL's, the relatively smaller effects including Zeeman cou- 
plings, disorders and Coulomb interactions are neglected. 
We obtain the flavor-dependent LL spectrum 
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where Uc = V^v/£b is the cyclotron frequency. The cor- 
responding eigenstates with n > can be formally writ- 
ten as a^^^^(t)n-i)^ for valley K and (/>n)^ 
for K'. For the anomalous n = LL's, the eigenstates 
are ((/)o,0)-^ with energy A and (0,^o)"^ with energy 
—A — 2Xsz' This clearly shows that the SU{4:) invariant 
four anomalous n = LL's are broken into a two- fold 
spin degenerate conduction band = LL at valley K 
and two spin split valence band n^^ = and n^^^ = LL's 
at K', as shown in Fig. [ijc), leading to integer quantum 
Hall effects at = and u = —1 but not u = 1. This 
is reminiscent to the anomalous n = LL's in few-layer 
graphene systems [28] . In graphene the SU (4) symmetry 
of n = LL's are completely lifted by electron-electron 
interactions [29ti3T| while the particle-hole symmetry re- 
mains, whereas in M0S2 trilayer both the SU{4:) and 
particle- hole symmetries are broken by the SOC and the 
inversion asymmetry. 

Other exotic features of LL's in M0S2 trilayers can be 
visualized in Fig.jljc) and further understood by expand- 
ing Eq. ^ Sit nB < 50 T: En,a = 2XrzSzSa- + a{A + 
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FIG. 1. (color online), (a) and (b) Electronic band structure 
near valley K and K'. (c) LL's with n = 0, 10, 80 orbitals. 
(d)-(f) Enlarged view of the LL's in Group I, II, and III in (c). 
The n 7^ LL's are broken into TzSz — 1 doublets and TzSz — 
— 1 doublets, due to the SOC and inversion asymmetry. LL 
crossing occurs between group II and III. The — LL is 
spin degenerate and only appears at valley K. The rijj , rijjj = 
LL's are spin- filtered and appear only at valley K'. 



A-Ar 8 with a = ±. (i) Because of the heavily 

massive Dirac Fermion character, the LL energy depen- 
dence on B appears to be linear rather than \fB. (ii) 
SOC break the LL's into two T invariant groups, with 
TzSz = 1 and TzSz = —1. However, each n 7^ LL is still 
doubly degenerate in each group, consisting of one spin '\ 
state from one valley and one spin ^ state from the other 
valley, (iii) The energies of two group LL's in the valence 
band not only have different slopes in B but also shift 
rigidly at 5 = 0, leading to LL crossing effects [32 at 
magnetic fields that satisfy 



Br = 



4A(A + A) 



(4) 



where n^j and n^^^ are the LL orbitals for Group II (r^s^ = 
1) and Group III {tzSz = —1) shown in Fig. [TJc). The 
second term is negligible when n^^ > 2.7n^^^. 

As an example, in Fig.|2]we consider the crossings be- 
tween LL's with = 38,39,40,41 from group II and 
LL's with njjj = 0, 1, 2 from group III. The n^^^ = LL is 
non-degenerate while other LL's are all doubly degener- 
ate. As a consequence, each region bounded by three 
LL's above the n^^j = LL has an odd filling factor 
v = — 2(njj + njjj) — 1 while bounded by four LL's below 
has an even filling factor u = — 2(njj +njjj, where n^j and 
n^jj are the orbital indices of the right and the lower LL's. 
These filling factors are all negative because the crossings 
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FIG. 2. (color online). Valence band LL crossing effects be- 
tween LL's with rijj = 38, 39, 40, 41 from group II (black) and 
LL's with rijjj = 0, 1, 2 from group III (red). All LL's are dou- 
bly degenerate except that the rZjjj = LL is non-degenerate. 
The negative number denotes the filling factor of the region 
bounded by the surrounding LL's. The crossing points have 
degeneracy p = 3 at rijjj = and p = 4 otherwise. 



only occur in the valence band. In addition, the cross- 
ing points at the n^^^ = LL all have degeneracy ^ = 3, 
while other crossing points all have degeneracy g = 4. 
The crossing of two LL's results in increased degenera- 
cies and lead to pronounced peaks in the measurement of 
longitudinal magnetoresistance, whose height is propor- 
tional to its degeneracy. Similar LL crossing effects in 
ABA-stacked trilayer graphene have been observed [32] 
via Shubnikov-de Haas oscillations. 

Even in the absence of interactions, the unconventional 
Hall plateaus follow the integer sequence u = ■ ■ ■ ^ —2M — 
6, -2M - 4, -2M - 2, -2M - 1, • • • , -5, -3, -1, 0, 
2, 4, • • • , where M = [4A(A + X)/{eBnv'^)] [33] reflects 
the fact that the n^^^ = LL lies between the LL's with 
rijj = M and n^^ = M + 1. The step of two in the 
sequence is a consequence of the r^s^ = ±1 classification, 
a hallmark of SOC. A step-one jump reflects the filling 
of a n = LL. The presence of u = and —1 arises 
from inversion asymmetry that makes the Dirac fermions 
massive and from SOC that rigidly shifts the two groups 
of valence band LL's with different TzSz characters. The 
switching between even and odd filling factors is a direct 
result of broken SU (4) symmetry among the anomalous 
n = LL's. We address that the n^i^n^u — LL's are 
spin-filtered and the exotic LL sequence is purely due to 
the non-interacting band structure. 

Spin-resolved p-n junctions. — A M0S2 P-n junction 
can be realized by using electrostatic gating to locally 
control the carrier type and density in two adjacent re- 
gions. In such a device, transport measurements in the 
quantum Hall regime reveal new plateaus with integer 
and fractional filling factors of two-terminal conductance 
across the junction. This effect arises from the redis- 
tribution of quantum Hall current among spin-resolved 



edge channels propagating along and across the junction. 
When Sz is a good quantum number, because the edge 
channels of n^i^n^^^ = LL's are spin-filtered, the full 
equilibrium must be achieved within each spin species 
separately. Consequently, the net conductance (in units 
of e^ /h) across the junction is quantized as follows 

Gpp^nn = min{|i/i^|, \u2t\} + min{|i/i;|, |i^2i|} , 

_ |^ltll^2tl , I^UII^2;| . . 

where + = vi and i^2t + ^2^ = ^2- In this limit, 
for a junction with vi = 2n and V2 = —1, electrons with 
one of the two spin-flavors are forbidden to flow into the 
p-doped region, since the only available LL n^j = is 
spin-filtered. Therefore, the net conductance is given by 
Gpn = n/(n + 1). The conductance across the junction 
becomes spin-independent when Sz is no longer a good 
quantum number, e.g.^ due to strong magnetic disorders. 
This limit is similar to the case of graphene [34 - 36 where 
all possible filling factors are even numbers because of 
the spin degeneracy. Therefore, the corresponding con- 
ductance reads Gpn = |^i I |^2|/(|^i | + |^2|) in the bipo- 
lar regime or Gpp^nn = niin{|i^i|, |i^2|} in the unipolar 
regime. Take the same example in which ui = 2n and 
V2 = —1, the net conductance in this spinless limit be- 
comes G = 2n/{2n + 1) instead. 

In addition to the novel transport properties, STM 
probes at the interface of two junctions can also detect 
a special fingerprint of the spin-filtered n^j = LL, re- 
vealing the strongly broken symmetry between the two 
groups of LL's in the valence band. As shown in Fig. [3j 
when the Fermi energy of one region is fixed between the 
two valence band n = LL's (region C), namely, between 
2A- A and -2A- A indicated by Eq. (|3|, while the Fermi 
energy of the other is outside this energy window, there 
will be one spin-filtered chiral edge state, among all the 
1^1 — ^2! channels, propagating along the interface. The 
chiral current will be controllable in the following senses, 
(i) Switching the magnetic field direction flips the spin- 
polarization of the current, (ii) Interchanging vi and 1^2 
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FIG. 3. (color online). Left panel: The solid (dashed) curves 
represent spin t (i) bands and the parallel lines denote their 
LL's. Regions A-D are energy windows separated by the three 
lifted energies of the four n = LL's that are depicted by red 
lines. Right panel: a schematic p-n junction with two different 
filling factors vi and 1^2 at two regions. 



switches the current direction while tuning ui and 1^2 ad- 
justs the current amphtude. (iii) Switching one of the 
Fermi level between A/B and D regions while fixing the 
other at C region changes the carrier type and fiips the 
spin-polarization. 

Spin Hall conductivity. — Clearly shown in FigjS) the 
spin Hall conductivity is quantized to asH = /h when 
the Fermi energy lies in the energy window C, due to 
the filling of spin-filtered n^j = LL. We address that 
(^SH does not vanish even in the absence of fields [37] - 
[40] , which is another consequence of the nontrivial va- 
lence band structure of M0S2 trilayers. In a massive 
Dirac fermion model, the momentum space Berry curva- 
ture [271 HI] ill the valence band is nontrivial and reads 

where m = A — XtzSz is the fiavor-dependent mass. 
At zero temperature, we obtain the valence band spin 
Hall conductivity (in units of e'^/h) by integrating 
Sz^zik^Tz, Sz) over the occupied states and summing 
over the spin- valley fiavors. As it happens in M0S2 
trilayers, the inversion asymmetry dominates the SOC 
(A > A) and the system is pinned to a QVH insula- 
tor in which asH = when the Fermi energy Cf lies in 
the gap, implied by the dependence of Qz- Due to 
the mass difference between TzSz = 1 and TzSz = — 1 
groups, however, ctsh contributions from the two groups 
are not completely canceled out when e/ crosses valence 
bands. These features can be seen in the u = line 
trace in Fig. [4]^a). If in the opposite limit, assuming 
A = 0, the system becomes a quantized spin Hall (QSH) 
insulator with asH = 2 as shown by the middle zone in 
Fig.[4]^b). Therefore, if the inversion asymmetry could be 
compensated externally, asH will be enhanced as long as 
Cf crosses the valence band, even though the system re- 
mains a QVH insulator. The inversion asymmetry arises 
from the difference between dz^ and d^id orbitals, which 
can be possibly modulated by applying an electric field 
or chemical doping, straining or electric gating the two 
S layers. To investigate the tunability of asH by varying 
Cf and by reducing inversion asymmetry, we propose a 
linearized phenomenological theory, i.e., the momentum- 
independent terms in Eq. ([T]) is replaced by 

u 

V{tz, Sz) = A(l - -^)rzSz + (A - u)az - XtzSzCTz , (7) 

where u is the modeled external potential that reduces 
the inversion asymmetry between dz^ and d^id orbitals. 
In the limit of u = this model recovers Eq. ([T]). In the 
special case of = A, inversion symmetry is restored and 
both az and TzSz terms vanish. 

Fig. [4] plots asH as a function of ej and u. The yellow 
(orange) line maps out where e/ touches the top of the 
valence band with r^s^ = +1 (r^s^ = — 1). In region 
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FIG. 4. (color online), (a) Valence band spin Hall con- 
ductivity as a function of e/ and u that reduces inversion 
asymmetry. The yellow (orange) line denotes the top of the 
TzSz = +1 ( — 1) valence band, (b) Quantum phase transitions 
between QVH and QSH states at = A ± A. 

C where e/ is inside the bulk gap, ctsh is identically 0. 
When Cf lies between the two valence band tops (region 
B), cTsH increases with decreasing ej. After Cf drops far 
below the lower valence band top (region A), asH starts 
to decrease slowly. Overall, asH is maximized when e/ is 
near the lower valence band top and grows substantially 
as u increases. When u is further increased, continuous 
quantum phase transitions (QVH-QSH-QVH) occur at 

= A^ A, as shown in Fig. [4]^b). Near the first (second) 
critical point, the mass of the Dirac fermions with TzSz = 
1 {^z^z = —1) changes sign and asH jumps from to 2 
(2 to 0) when e/ is inside the gap. 

Discussions. — The Hall plateau in graphene follows 
the sequence v = 4(n + 1/2), with SU{4:) symmetry 
breaking only in the case of high fields, weak disorders, 
and strong interactions |29H3T] . In a M0S2 trilayer, quan- 
tum Hall ferromagnetism arises naturally even in the ab- 
sence of interactions. SOC and inversion asymmetry also 
imprint a single-particle signature on the LL spectrum: 
the broken symmetry of the n = LL's and their en- 
ergies are independent of the field strength. The un- 
conventional Hall plateau sequence in a M0S2 trilayer 
becomes richer in a p-n junction with the appearance 
of a spin-resolved fractionally quantized conductance. 
A controllable spin-polarized current materializes within 
this geometry and can be probed by STM at the inter- 
face. Consequently, the MoS2 trilayer may find use in a 
multi-functional spintronic device. Like in ABA trilayer 
graphene [32 , LL crossing effects occur in the M0S2 tri- 
layer but only in its unconventional valence band. The 
two valleys are imbalanced in their contributions to ctsh^ 
providing hope for increasing asH by reducing the gap. 
With improvements on its mobility, the M0S2 trilayer 
system may realize as an alternative to graphene in ful- 
filling the desire for a gapped Dirac system with strong 
SOC, and we also anticipate observations of its unique 
Hall phenomena discovered in this Letter. 
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